一维流体饱和粘弹性多孔介质层的动力响应

本文研究了不可压流体饱和粘弹性多孔介质层的一维动力响应问题。基于粘弹性理论和多孔介质理论，在流相和固相微观不可压、固相骨架服从粘弹性积分型本构关系和小变形的假定下，建立了不可压流体饱和粘弹性多孔介质层一维动力响应的数学模型，利用Laplace变换，求得了原初边值问题在变换空间中的解析解，并利用Laplace逆变换的Crump数值反演方法，得到原动力响应问题的数值解。数值研究了饱和标准线性粘弹性多孔介质层的动力响应，分析了固相位移、渗流速度、孔隙压力及固相有效应力等的响应特征。结果表明，与不可压流体饱和弹性多孔介质相同，不可压流体饱和粘弹性多孔介质中亦只存在一个纵波，并且固相骨架的粘性对动力行为有显著的影响。     

泄爆过程中二次爆炸的动力学机理研究

在容积为0.00814m3的柱形泄爆容器中，对泄爆现象进行实验研究. 容器内充满当量比为1的甲烷-空气预混气，采用底端中心点火，泄爆压力为230±15kPa. 基于k-ε湍流模型和EBU燃烧模型，利用同位网格的SIMPLE算法，对该现象进行了数值模拟. 实验和计算获得的外轴线上4个测压点的压力曲线和外流场的阴影和数值照片，形象地描述了高压泄爆时外部流场的变化. 数值结果与实验结果基本一致. 根据实验和数值结果，详细地讨论了泄爆过程中二次爆炸产生的动力学机理. 泄爆的初始阶段，在破膜激波的引导下，泄出的未燃气体因欠膨胀在外流场形成稀疏波低压区和悬激波高压区. 高压区可燃气体密度和温度上升，成为高密度的预热区域. 随后，火焰以射流形式从泄爆口泄出，点燃可燃气云. 受湍流等因素的影响，特别在高密度的预热区域，燃烧速率可能迅速增大，从而导致二次爆炸.     

一类生物流体力学模型的稳定性

在研究颅内压集中参数模型的基础上，改进自适应的龙格-库塔法对一类生物流体力学模型进行数值模拟。通过合理控制计算量，得到了微分方程近似解的局部截断误差的估计。使用矩阵特征值分析微分方程的稳定性，在实际生理范围内变化模型参数，得到特征值随参数变化的关系，指出模型中存在实部为正的特征值。文章还讨论了控制矩阵特征值的变化对数值方法稳定步长的影响，并得到步长的取值范围。通过理论分析。得到控制方程的解随时间的发展和模型中生理参数的选取相关。分析了特征值变化的稳定性和参数的关系，并对模型进行改进，讨论其稳定性的情况。     

Direct Numerical Simulation of Statistically Stationary One- and Two-Phase Turbulent Combustion: A Turbulent Injection Procedure

A new turbulent injection procedure dedicated to fully compressible direct numerical simulation (DNS) or large eddy simulation (LES) solvers is proposed. To avoid the appearance of spurious acoustic waves, this method is based on an accurate tracking of the turbulent structures crossing the boundary at the inlet of the domain. A finite difference DNS solver has been coupled with a spectral simulation in which a statistically stationary homogeneous turbulence evolves to provide fluctuating boundary conditions.A new turbulence forcing method, dedicated to spectral solvers, has been developed as well to control the major properties of the injected flow (turbulent kinetic energy, dissipation rate and integral length scale). One-dimensional Navier–Stokes characteristic boundary conditions extended to non-stationary flows are coupled with the injection procedure to evaluate is potential in four various configurations: spatially decaying turbulence, dispersion of vaporizing sprays, propagation of one- and two-phase V-shape turbulent flames.     

水中爆炸冲击波作用下混凝土墩动态响应初步分析

用水下爆炸冲击分析（ＵＳＡ）方法,在假设条件下,数值模拟了混凝土结构物水中爆炸冲击波作用下的加速度及压力响应,是对浅层水中爆炸流－固耦合动力学问题数值解的初步尝试。     

Effect of the Geometry of the Working Chamber of a Cylindrical-Conical Hydrocyclone on the Non-Newtonian Fluid Dynamics

Numerical simulation is used to study the dynamics of non-Newtonian free-surface flow in a cylindrical-conical hydrocyclone. For different angles of taper of the conical section of the hydrocyclone, the pressure and velocity distributions are calculated, together with the dependence of the fluid film thickness on the axial coordinate. The effect of the rheological properties of the fluid and the controlling similarity parameters on the flow dynamics is studied.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, 2005, pp. 102–112.Original Russian Text Copyright © 2005 by Yablonskii.     

Averaging methods for numerical simulations of flows through heterogeneous porous media

Many natural rock systems contain small patches of different permeability which affect the flow of fluids through them. As these heterogeneities become smaller and more numerous, they become harder to model numerically. We consider how to reduce the computational effort required in simulations by incorporating their effects in the boundary conditions at the edges of each grid block. This is in contrast with current methods which involve often arbitrary changes in the fluid properties. The method is restricted to the case of widely-spaced patches, which simplifies interaction effects. The system then reduces to an array of dipoles, and two averaging methods are proposed for finite grid blocks. Several infinite systems, including vertical and horizontal bands, are also considered as further approximations. There is a great wealth of existing results from different fields which lead to identical mathematical problems and which can be used in these cases. Finally, we consider how to use these techniques when the precise configuration of the grid block is not known, but only its statistical properties. This can lead to results which are very different from the deterministic case.     

用HLLC方法处理运动边界

针对含有运动边界非结构网格点的显式有限体积方法，研究了如何使用HLLC方法处理运动边界条件，计算了活塞拉伸、压缩问题，并和解析解进行了比较．作为应用示例，比较了静止圆柱绕流和静止流场中运动圆柱的绕流情形。结果表明：应用HLLC方法处理运动边界是可行的。     

Effect of Boundary-Layer Thickness on the Structure of a Near-Wall Flow with a Two-Dimensional Obstacle

Results of physical and numerical experiments on investigating the effect of the depth of immersion of a two-dimensional obstacle with a square cross section into a developed turbulent boundary layer on the length of the separated flow region are presented. The numerical simulation is based on solving averaged Navier–Stokes equations with the use of the k– model of turbulence. The near-wall flow is visualized in the experiments, and the fields of mean and fluctuating velocities are measured. Flow regions where the results of numerical simulation agree with experimental data are determined. It is shown that the length of the recirculation flow region in the near wake increases with decreasing depth of immersion of the two-dimensional obstacle into the turbulent boundary layer.     

On numerical solution of dynamic problems in the theory of reinforced shells

Dynamic problems for cylindrical shells reinforced with discrete ribs are examined. A numerical algorithm based on Richardson extrapolation is developed. Specific problems are solved, and the results are analyzed __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 5, pp. 50–56, May 2006.